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ABSTRACT 


Polarization, differential cross-section, and total reaction cross-section data for 
the scattering of 16- to 22-MeV protons from complex nuclei were analyzed with the op- 
tical model. An average potential was obtained which describes most of the data. De- 
viations from the average spin-orbit potential were observed for nuclei with doubly 

1 fi 

closed shells. The analysis of the energy dependence of proton scattering from O 
showed that a smooth variation of the optical model parameters results in satisfactory 
fits to the data between 30 and 50 MeV. The parameters cannot be extrapolated into the 
20-MeV range. 



OPTICAL MODEL ANALYSIS OF PROTON SCATTERING 
IN THE RANGE 16 TO 22 MeV 


by Edmund T. Bosch Itz 
Lewis Research Center 

SUMMARY 

An optical model analysis was performed of differential cross sections, polariza- 
tions, and total reaction cross sections for the scattering of protons from many nuclei 
in the energy region 16 to 22 MeV. An average potential consisting of nine parameters 
was obtained which predicts polarizations, differential cross sections, and reaction 
cross sections fairly well for most nuclei with masses larger than 30. 

The parameters of the real central potential agree with those proposed earlier by 
Perey. The strength and width of the surface absorption potential increase with the tar- 
get mass. Fits to the polarization data yield a spin-orbit potential with a radius and dif- 
fuseness smaller than those of the real central well. Anomalously low values for the 

parameters of the spin-orbit term were found in the mass regions of doubly closed shells. 

40 

The effect is most pronounced for Ca . The analysis of the energy dependence of the 

1 fi 

proton scattering from O showed that it is impossible to extrapolate the smooth varia- 
tion of the optical model parameters between 30 and 50 MeV to lower energies. 


INTRODUCTION 

It is an established fact that the elastic scattering of intermediate energy nucleons 
by heavier nuclei can be well described by the nuclear optical model. Most of the analy- 
ses, however, until recently, were based on differential cross-section data only, since 
polarizations and total reaction cross sections were rare. For this reason, the analyses 
gave little information about the spin -dependent part in the potential and left uncertainties 
in the imaginary part. 

Systematic polarization studies were started by Rosen and Brolley (ref. 1) in 1958. 
Their extensive work between 8. 0 and 14. 5 MeV covered a wide range of nuclei. Unfor- 


tunately, the polarization observed at these energies becomes quite small for nuclei in 
the upper half of the periodic table. 

The elastic scattering of 18. 5 and 20. 5 MeV polarized protons by complex nuclei was 
studied extensively at this laboratory in past years (refs. 2 to 4). The experimental 
work was motivated by the need for more systematic polarization data at higher energies 
where sizable polarizations were expected for the heavy nuclei as well. In addition, the 
scattering by sets of isotopes and isobars should reveal the sensitivity of the polarization 
to the nuclear surface structure. This is quite plausible because the spin-orbit potential 
in its presently accepted form is proportional to the derivative of the central potential. 
Light nuclei were included in the polarization experiment in order to learn about the lim - 
itations of the optical model. 

Parallel to this work, Baugh et al. (refs. 5 and 6) have investigated the polarization 
for 17. 8 -MeV protons scattered from a number of nuclei. The measurements extend 
only to a 120° laboratory angle. This is unfortunate since the polarization distribution 
at large angles is quite significant for the analysis. The results of an optical analysis 
of data limited in angle might contain additional ambiguities. 

More recently polarization data have become available at 30 and 50 MeV from the 
Rutherford Laboratory (refs. 7 to 10) and UCLA (private communication from W. T. H. 
van Oers), and at 40 MeV from the University of Minnesota (ref. 11) and the Oak Ridge 
National Laboratory (ref. 12). 

In the present report, an optical model analysis is reported which uses polarization 
and differential cross-section data from this laboratory and also cross sections from the 
literature. In this study the emphasis is placed on the properties of the spin-orbit po- 
tential (i. e. , its shape and its dependence on the bombarding energy and on the mass of 
the target nuclei). 


SYMBOLS 


A nuclear mass number 

a diffuseness parameter of real central potential 
aj diffuseness parameter of surface imaginary potential 
a gQ diffuseness parameter of spin-orbit potential 
c velocity of light 

E proton energy 

F unit of length (10 - ^ cm) 

h Planck’s constant divided by 2n 


2 


E-4583 



I 

m 

77 

P 

exp 
AP 


exp 


th 


r 

r 


so 


V 

V c 

V so 

W 

W 

Z 

0 

CT 

CT 

exp 
Act 


D 


CM 


exp 


th 


angular momentum vector 

rest mass of tt meson 

experimental polarization 

estimated error in P gX p 

theoretical polarization 

radius of real central potential 

radius of Coulomb potential 

radius of imaginary potential 

radius of spin -orbit potential 

real central potential 

Coulomb potential 

spin -orbit potential 

volume part of imaginary potential 

surface part of imaginary potential 

nuclear charge number 

scattering angle in center of mass system 

Pauli spin operator 

experimental differential cross section 
estimated error in a eX p 
theoretical differential cross section 

measure of goodness of fit 


PREVIOUS ANALYSES IN THE RANGE 16 TO 22 MeV 

Two detailed optical model analyses exist in the literature in the energy range of 
the NASA data. Perey (ref. 13) has analyzed the differential cross sections for elastic 
protons scattering between 7 and 22 MeV with the emphasis on finding smooth trends in 
the parameters. His analysis resulted in an average potential which contains energy and 
mass dependences in a number of parameters. Rosen et al. (ref. 14) on the other hand, 
used the large amount of polarization data at 8. 0, 10. 5, and 14. 5 MeV in addition to the 
existing cross-section data in order to establish a simple M universal M potential which 
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TABLE I. - COMPARISON OF "AVERAGE" POTENTIALS GIVEN 


BY PEREY (REF. 13) AND ROSEN ET AL. (REF. 14) 


Parameter 

Rosen 

Perey 

V, MeV 

53.8 - 0.33E 

53.3 - 0. 55E + 0.4 Z + 27 N ' Z 



A 1/3 A 

W, MeV 

7. 5 

13. 5 ±2 below E p = 17; 3A 1 / 3 ± 1. 5 above 

V so , MeV 

5. 5 

8. 5 for 17 — Ep — 22; 7. 5 for E p < 17 

a = a so> F 

0.65 

0.65 

®T F 

0.70 

0.47 

r = r T = r cn , F 
o I so* 

1.25 

1.25 


would describe the scattering of nucleons by nuclei over a wide range of energies and 
nuclear masses. 

The two sets of parameters are listed in table I. The various terms are defined in 
the symbol list. The simplicity of Rosen’s potential is quite appealing. In figure 1 pre- 
dictions from both potentials are compared with the measured polarizations at 20. 5 MeV. 
The solid lines correspond to Rosen’s potential, and the dashed ones to Perey's para- 
meters. Both potentials do poorly on the light nuclei. For the medium weight nuclei, 
Rosen’s predictions are fairly good, but with increasing target mass the theoretical pre- 
dictions become too large and get ’’out of phase. ” Similar difficulties are encountered 
in fitting the differential cross sections at 22. 2 MeV. Perey’s potential shows the oppo- 
site trend; it does poorly in the nickel region but improves for the heavier nuclei. The 
second fact can be attributed to the increase of the imaginary well depth with target 
mass. 

In regard to predicting the total reaction cross sections (refs. 15 and 16), Rosen’s 
potential produces a value which is 30 percent too large for nuclei up to mass 100; 

Perey’ s potential gives reaction cross sections which are too low for heavy nuclei (see 
fig. 6). 
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Rosen's analysis 

Perey's analysis 

• Measurements made by 

Bercaw et al. (ref. 20) at 20 MeV 



Center of mass angle, 0 q M , deg 


Figure 1. - Polarizations predicted by potentials of Rosen and Perey. 




METHODS OF ANALYSIS 


The present analysis was performed with a modified version of SCAT 4 (ref. 17) and 

o 

Davidon’ s automatic search routine (ref. 18), which minimizes the quantities X and 

2 CT 
Xp where 


n 


X 2 = 

a 


x 2 - 

A p " 


E 

i=l 

n 

E 


i=l 


Mi> - CT exp( 0 i)' 


n 2 


Ao exp^P 


FW - P exp( 0 i) 


,1 2 


P exp^P 


The calculation was performed on the IBM 7094 computer. Although the minimum 

2 

X was generally accepted as a criterion of best fit, visual inspection of goodness of 

2 

fit served as an additional requirement. This becomes necessary when X D is large 

2 v 
and a decrease in XT, contradicts the shape of the polarization angular distribution. 

r 2 

Also, an inferior fit from the X point of view might be acceptable when the angular 
resolution typical for polarization measurements is considered. The optical potential 
used in this analysis is local and has the usual form 


V(r) = V c 



i|w-4W D 





l 


where 


x = 


r - r A 
o 


1/3 


r 

x T = - 


ryP / 3 

a I 


X 


so 


r - r A 
so 


so 


1/3 
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and V is the Coulomb potential corresponding to a uniformly charged sphere of radius 
r c = 1.25 A 1//3 F. 

At the time when the polarizations were measured at this laboratory, it was assumed 
that later differential cross sections could also be measured for all nuclei at the same 
proton energy. Unfortunately - until now - matching cross sections became available 
only in a few cases. In most cases one had to be content with differential cross sections 
and total reaction cross sections at neighboring energies and in some cases even at 
neighboring nuclear masses. The following data were used: 

(1) 40 polarization angular distributions (between 15° and 165° laboratory angles) of 

the energies 16. 6, 18. 5, and 20. 5 MeV (Bercaw et al. (refs. 19 and 20)) 

(2) 25 differential cross sections at 17. 0 (Dayton and Schrank (ref. 21)), 18. 5 (Eccles 

et al. (ref. 22)), 20. 5 (Gray et al. (refs. 23 to 25) and private communication 

from N. Baron, Lewis Research Center), and 22.2 MeV (Fulmer (ref. 26)) 

(3) 7 total reaction cross sections at 16. 4 (Pollock and Schrank (ref. 15)) and 

22. 2 MeV (Fulmer (ref. 16)) 

The analysis of such unrelated sets of data required a special approach. The analy- 
sis was performed in two steps: 

(1) In the first step an attempt was made to find an improved "average" potential 

(with constrained spin -orbit coupling, i. e. , a gQ = a and r g0 = r Q ) which would give a 

"reasonable" fit to all data in this energy range. This was done in various ways. In 

faborable cases where polarizations and differential cross sections existed at the same 

energy, the sum XZ, = X~ + XZ, was minimized directly. In other cases where both 

O' sr 2 2 2 

kinds of data were available only at different energies, either X~ = X + Xt, was mini- 

mized for the average energy or X* and Xp were minimized separately (each at the 
correct energy); then the "optimum" parameter pairs were gradually adjusted to con- 
verge to one value which produced a compromise fit to both kinds of data. In the major - 

2 2 

ity of the cases, the best method appeared to be the minimizing of X^. and X p for all 
nuclei and energies individually for a given search on two or three parameters. The 
widely fluctuating parameters which resulted from this search were then plotted against 
nuclear mass and bombarding energy and common trends were investigated. These 
trends were used in the next search on different parameters and new smooth trends were 
established. This method seems appropriate since one is looking for "average" trends 
in the true spirit of the optical model. There are much less data available for the total 
reaction cross sections. Therefore, they were not included in the search, but they 
served as a guide to reject parameter values which did not predict the reaction cross 
sections correctly. 


7 


(2) In the second step, the trends of the spin-orbit potential were to be investigated. 
In order to see those trends, the fluctuations of the other parameters must be averaged 
out. This was done by using the average central potential from step (1); now the search 
routine varied V go , a gQ , and r gQ , minimizing Xp only. An analysis under such "re- 
stricted conditions" is considered permissible since it is essentially the polarization 
which defines the parameter space of the spin-orbit potential. This assumption has been 
confirmed at low energy by Goldfarb et al. (ref. 27). They have shown that small devia- 
tions in the spin-orbit parameters could be treated as perturbations which affect the 
polarization but hardly change the cross section. 


SEARCH FOR AVERAGE POTENTIAL 
Central Potential (Constrained Spin-Orbit Potential) 

The search was started from Rosen’s potential since it is based on the analysis of 
both cross-section and polarization data. Also the simplicity of this potential is appeal- 
ing. With the geometrical parameters kept fixed at first, combined searching and rough 
gridding on the depths of the potentials V, W D , and V go (W was set equal to zero) 
showed the following trends: All searches on W D clearly indicate the need for an in- 
crease with nuclear mass. This mass dependence reduces the magnitude of the polari- 
zation for heavy nuclei, which is the main drawback of Rosen' s parameter set, but it 
does not allow for the difference in the "phase” between theoretical and experimental 
polarization distributions. This difference can be removed only by a mass dependence 
in the real central well. Studying the mass (and energy) dependence of V in detail re- 
vealed that both dependences were reasonably close to the ones given by Perey (ref. 13). 
Therefore, from that point on the real potential V was set equal to Vp ere y- Searches 
on V go (a gQ = a, r gQ = r Q ) confirmed Rosen’s value of 5. 5 MeV. Now V, V go , a = a gQ , 
r 0 = r go were kept fixed, and the imaginary potential was inspected more closely. In 
this part of the analysis, the total reaction cross section served as "a guide" to exclude 
sets of parameters which wbuld give good fits to polarization and differential cross 
section, but which would not agree with the total reaction cross sections. In order to 
obtain the correct mass dependence of the imaginary well, the radius parameter rj was 
varied in steps of 0. 05 F between 1. 15 and 1. 35 F. For each value of rj, the diffuse- 
ness aj was varied in steps of 0. 05 F between 0. 45 and 0. 75 F. For these various 

combinations, the computer code searched on W n . The parameters r T = 1.25 F, 

1/3 u 1 

aj = 0. 65, and W^ = 2A ' were found as an optimum set. Yet, the prediction of the 

reaction cross sections still needed improving. In reference 28 it is suggested that a.^ 

might be larger for heavy nuclei, so various mass dependences were tried for aj; 
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aj = 0. 54 + (A/1000), aj = 0. 45 + (A/250), aj = 0. 38 + (A/550), and aj = 0. 45 + (A/650) 
were chosen. The last mass dependence is most successful, and the corresponding 
values for W D indicated that W = 1. 5A 1//3 + 2. 5. The inclusions of a volume absorp- 
tion W had little effect on the reaction cross section; it increases only by about 1 per- 
cent per 1 MeV well depth. In presenting average parameters for the imaginary poten- 
tial, it should be pointed out that there are large fluctuations in W^, as shown in fig- 
ure 2. These fluctuations are probably due to the fact that contains the " individu- 

alism” of the nucleus in regard to nuclear reactions and modes of excitation. 
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Figure 2. - Optimum values of W D obtained from searching on 
Wq with remaining parameters fixed, (aj = 0.65 F.) 


Independent Spin-Orbit Potential (Fixed Central Potential) 

After an average central potential was found, V, a, r Q , W D , aj, and rj were fixed 
and the search code varied the parameters of the spin-orbit part, minimizing Xp only. 
The optimum values from this "restricted” analysis were checked in eight cases where 
matching cross-section data were available. In these cases the search code varied V 
and W in addition to V go , a gQ , and r go , which now minimized the sum X 3 , = Xp + X 3 . 
The optimum values for the spin-orbit potential parameters found by this "unrestricted" 
search differed only insignificantly from those obtained under the restricted condition, 
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verifying the validity of the method. The values for V__, a 0( _, and r__ which corre- 
spond to the smallest Xp are plotted as a function of nuclear mass in figure 3. This 
plot contains the results of the analyses of both 18. 5- and 20. 5-MeV data, since no 
energy dependence of the average parameters was apparent. With the exception of the 
light nuclei, the average trend of the parameters is indicated by the horizontal lines with 
the range of scatter given by the crosshatching. The average values are V go = 5. 5±0. 5 
MeV, a go = 0. 55±0. 05 F, and r gQ = 1. 12±0. 05 F. 


Average trend of parameters 

Range of scatter 

• • 



A 


Figure 3. - Optimum values of V s0 , a S p, and r so obtained 
from searches on these parameters with parameters of central 
potential fixed. (V = V Pprpv , W - 2A^, a = a T = 0. 65 F, and 
r = r r 1.25 F.) V 

These average values together with the average parameters for the central potential 
(V = Vp e , a = 0. 65 F, r = 1. 25 F, W D = 1. 5 A 1 / 3 + 2. 5, aj = 0. 45 + (A/650), and 
rj = 1. 25 F) produce the theoretical curves in figures 4 to 6. The overall agreement 
with the experimental data is quite good. The individual fits can certainly be improved if 
the well depths V, W^, and V go can vary freely. For a quantitative comparison, the 
paremeter sets of Rosen and Perey are listed in table II together with the average para- 
meters from this analysis. Considering recent polarization data, Perey modified his 
first potential (ref. 29). The revised parameters are listed in column 4. In columns 
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Polarization, 


I 




(a) Proton energy E p = 18.5MeV. 




Figure 4. - Polarizations predicted by "average” potential listed in table II, column 6. (Data are from refs. 2 and 20. ) 


li 



5 and 6 the results from this analysis are presented. The set in column 5 produces good 
fits to the polarizations and differential cross sections but is inferior in predicting the 
total reaction cross sections. The set in column 6 corrects for this. The comparison 
shows fairly good agreement with Perey’s analysis, the difference being in Wj^, aj, 

V go , and a g0 . The mass dependence of V, W D , and aj are essential features of the 
optical model. Neglecting this dependence in Rosen’s potential causes inferior fits to 
the data. 






(b) Proton energy E p = 20. 5 MeV. 
Figured -Concluded. 
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Differential cross section, do/dQ, mb/sr 





TABLE H. - COMPARISON OF PREVIOUS POTENTIALS WITH 


BEST "AVERAGE" POTENTIALS FOUND IN THIS ANALYSIS 


Parameter 

Rosen 

Perey 

Present work 



First 

Second 

First 

Second 

V, MeV 

53.8 - 0.33E 

q _ n i n a ^ . on N - Z 


IT 

IT 

A l /3 a 

V P 

V P 

V P 

W D , MeV 

7. 5 

3A X /3 

3A 1 / 3 

CO 

> 

a 

1. 5A 1 / 3 + 2. 5 

V Qrt , MeV 
so’ 

5. 5 

8. 5 

6.0 

5. 5 

5. 5 

a, F 

0.65 

0.65 

0.65 

0.65 

0.65 

a so> F 

0.65 

0. 65 

0.47 

0.55 

0. 55 

a r , F 

0.70 

0.47 

0.47 

0.65 

0.45 + — 
650 

r o> F 

1.25 

1.25 

1.25 

1.25 

1.25 

r p F 

1.25 

1.25 

1.25 

1.25 

1.25 

r so> F 

1.25 

1.25 

1.12 

1. 12 

1. 12 

X 2 ( a 17. 0) 

a 9. 5 

6.4 

10.2 

6. 1 

8.0 

x| ( a 18. 5) 

27.0 

49.0 

15. 5 

11.6 

11.8 

X 2 ( a 22. 2) 

20.4 

16.7 

19.4 

6.3 

6.4 

x| ( a 20. 5) 

29.3 

32.0 

14.8 

5. 7 

6. 5 

X 2 

A av 

21. 5 

26.0 

15.0 

7.4 

8. 1 


Q 2 

a Values for X are number of X f s per data point. 

ANOMALIES IN THE SPIN-ORBIT POTENTIAL 
(THE NUCLEI K 39 , Ar 40 , AND Ca 40 ) 


In figure 3 the average trend of the parameters V gQ , a gQ , and r gQ is shown. 

There are deviations from the average values for the light nuclei (A < 30) which are dis- 
cussed later. There are other mass regions, however, where considerable deviations 
from the average values occur. Anomalously low values are found around A = 40 
(K 39 , Ar 40 , and Ca 40 ), A = 90 (Sr 88 and Zr 90 ), and A = 208 (Pb 206 and Pb 208 ) which 
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seem to be related to shell closure. The effect is most pronounced in the A = 40 region 
and less so in the other regions. The anomalous parameters show no regularities. For 
one nucleus, only r gQ shows irregularity; for another, V gQ and/or a gQ ; sometimes 
all three parameters have to be decreased to obtain good fits. This could be due to some 
ambiguities between V gQ , a gQ , and r gQ similar to the well-known ones between V and 
r and W D and a^. If some ambiguities between the parameters of the spin-orbit po- 
tential do exist, then it might be possible to keep one or two parameters fixed and to 
assign the anomalous behavior to the remaining parameters. This was tried, but gener- 
ally it was impossible to attribute the anomalies to just one or two of the three para- 
meters. Some progress, however, was made in the A = 40 region, which is described 
in the following section. 

Ca 40 has been an interesting nucleus for quite a while. In the contest of the optical 

model, it is interesting because of the difficulties encountered in fitting the experimental 

data In the course of the experiments at this laboratory, large differences in the polari- 

40 39 

zation distributions were observed in the scattering of 18. 5-MeV protons by Ar , K , 
and Ca 4 ® (ref. 30). 

Volkin performed an optical model analysis (ref. 31) of these data in order to explain 

these differences. He found that best fits were obtained by decreasing the radius para- 

40 39 40 

meter r from 1. 25 to 1. 18 F and 1. 03 F when proceeding from Ar to K and Ca . 

SO 

He concluded that this decrease might be indicative of a shell closure effect in the nu- 
clear optical model which appears in the form of the spin-orbit potential. For closed 
shell nuclei, the spin-orbit potential is centered approximately at the nuclear matter 
radius rather than at the central potential radius. Volkin also considered a complex 
spin -orbit potential to fit the data, but he concluded that there is no need for an imagi- 
nary part in the spin-orbit potential if a smaller radius (r__ < r ) is used. 

40 ® 

Baugh et al. (ref. 32) analyzed their p-Ca polarization angular distribution at 

17. 5 MeV with an imaginary part W ori = 2 MeV. To elaborate on this result, they meas- 
ured the energy dependence of the polarization for this nucleus between 14. 0 and 

18. 5 MeV (at four laboratory angles between 85° and 115°) and confirmed the earlier 

conclusion that W_. « 2 MeV over this energy range. 

6,0 39 40 

In the present work, the polarizations and differential cross sections for K , A , 

and Ca 40 were analyzed in the energy range 14. 5 to 21 MeV. At first, only V on , a. an , 

9 oU oU 

and r were varied in combined rough gridding and searching; XT, was thus minimized. 

SO v 

It was found, as mentioned previously, that all three parameters tend to have smaller 

2 

values than the average ones; but, if somewhat inferior fits (~30 percent increase in X ) 

2 2 2 

are accepted, r gQ = 1. 1 F can be kept fixed. In the following searches, X T = X p + X^ 
was minimized when V, W, V gQ , and a gQ were varied (where a = aj = 0. 65 F, r Q = 
rj= 1.25F, and r gQ = 1. 1 F). Both V and W were included in the search because 
differential cross-section data were used in this part of the analysis. The results are 
shown in figure 7; the corresponding parameters are listed in table m. A number of 
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interesting features should be pointed out. First, the depth of the real central well is 
40 40 

smaller for Ar than for Ca contrary to the predictions from the asymmetry term in 

Perey’s potential, where V . n > V 4n by 2. 7 MeV. No explanation is found for this 

Ar* u Ca 4U 

result. Second, the depths of the imaginary central well decrease monotonically with pro- 
gressive shell closure, Ar^® Ca^®. This might indicate that fewer inelastic channels 
are open for the closed shell nucleus. Third, the depths and widths of the spin-orbit 
potential, V go and a gQ , individually do not show a systematic behavior, but their pro- 
duct V go a go does. This quantity might be interpreted as the ” effective strength” of 
the spin -orbit coupling. The energy dependence of this ” effective strength” is plotted 
in figure 8. The crosshatched area marks the range of scatter of the average parame- 
ters. The interesting features of this energy dependence are the systematically lower 

40 

values of V_„a ort for Ca at all energies and the occurrence of a minimum in the spin- 
so so jq 4 Q 

orbit force for both Ar and Ca near 17 MeV. Both features might be understood 

40 

qualitatively. The lower values of V go a go for the Ca nucleus can probably be ex- 
plained in terms of the spin-dependent part of the nucleon -nucleon potential. The 
nucleon-nucleon potential contains a tensor term (of longer range) and a spin-orbit term 


TABLE m. - OPTIMUM PARAMETERS FOR Ar 40 , K 39 , AND Ca 40 OBTAINED 
FROM SEARCHES ON V, W D , V g0 , AND a g0 WITH THE USE OF 
POLARIZATION AND CROSS-SECTION DATA 


[Geometric parameters: a = a^ = 0. 65 F, r Q = rj = 1. 25 F, and r g(> = 1. 10 F.] 


Energy, 

Nuclei 

Depth of 

Depth of 

Depth of 

Width of 

Effective 

Goodness 

E, 


real 

volume 

spin- 

spin- 

strength, 

of fit, 

MeV 


poten- 

part of 

orbit 

orbit 

V a 
v so so 

2 

X 



tial, 

imagi- 

potential, 

potential, 


i 



V 

nary, 

v so 

so 






W 





14. 5 

Ar 

47.4 

7. 2 

5. 5 

0.55 

3.0 

159 


Ca 

48.4 

5.2 

5.0 

.37 

1.9 

412 

17.0 

Ar 

46.7 

7.7 

3.7 

0.43 

1.6 

550 


Ca 

47.9 

5.2 

4.2 

.25 

1.0 

361 

18.3 

Ar 

47.4 

8.5 

4.4 

0.60 

2.6 

357 


K 

48.7 

6.4 

4.0 

.49 

1.9 

744 


Ca 

48.0 

5.4 

3.8 

.34 

1.3 

724 

21. 0 

Ar 

43.5 

7.7 

4.8 

0.60 

2.9 

275 


K 

a 47. 1 

a 8.4 

4. 5 

.57 

2.6 

340 


Ca 

a 46. 2 

a 8. 2 

4.6 

. 53 

2.4 

249 


a Obtained from searches that used cross sections of neighboring nuclei. 
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Polarization, 



Center of mass angle, 0^, deg 
(3) Ar 40 . 

Figure 7. - Comparison of calculated polarizations and differential cross sections (parameters of table III) with experimental data of 
references 2, 20, 23 to 26, and private communication from N. Baron. 
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(b) K 39 . 

Figure 7. - Continued. 


40 

(of shorter range). For the doubly closed shell nucleus Ca , the tensor term will 
vanish, and the remaining (short range) nucleon -nucleon LS-term is responsible for the 
smaller parameter values in the spin -orbit term of the optical potential. The occurrence 
of a minimum in the ’’effective strength” near 17 MeV might be correlated with the en- 
tering of a new partial wave. In the scattering of 14. 5-MeV protons from complex 
nuclei, Rosen and Stewart (ref. 33) have observed that in mass regions where p . R = nK 
(corresponding to the entering of a new partial wave), the back angle polarization patterns 
are extremely sensitive to the addition of even a single nucleon in the target nucleus. 

The present analysis of A 40 , K 39 , and Ca 40 indicates a similar sensitivity of the polari- 
zation. The question arises if a decrease in the spin-orbit force for Ca 49 would also 
appear at the entering of the next partial wave at 30 MeV (p • R = 5fi). Recently, polar- 
ization and cross-section data at 30 MeV (Rutherford Lab. ) and at 40 MeV (Oak Ridge 
National Lab.) have been analyzed by Satchler (ref. 34) and Blumberg et al. (ref. 35). 

For Ca 49 considerable difficulties were encountered in obtaining satisfactory fits to 
polarizations and cross sections simultaneously. However, the best fits obtained at 
30 MeV required reduced values of V go and a g0 , resulting in a product of about 1 (MeV) 
(F), which is comparable to the findings at 17 MeV. At 40 MeV, no such effects are seen. 
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Figure 8. - "Effective strength" (V so a S g) of spin - 
orbit coupling as function of bombarding energy. 

LIGHT NUCLEI (ENERGY DEPENDENCE OF p-O 16 ) 

It has been shown that the average optical model potential (columns 4 and 5 of 
table II, p. 15) can describe differential cross sections, polarizations, and total reaction 
cross sections for the heavier nuclei very well. From figure 4 (p. 11) it is clear, how- 
ever, that this is not true for the light nuclei (A < 30). In some cases (C* 2 at 20. 5 MeV 
28 32 ~ 

and Si and S at 18. 5 MeV), there may be a qualitative agreement between theory and 
experiment; however, in other cases (N^ and O 1 ®), there is complete failure. If the 
optical model parameters are allowed to vary freely, good fits can generally be obtained, 
but the parameters vary rapidly with energy and from one nucleus to the other. There 
are certain trends in some of the parameters, as can be seen in figure 3 (p. 10). The 
parameters of the spin-orbit coupling, in particular, the diffuseness a , tend to smal- 
ler values, but the meaning of this is not clear. Other analysis in this energy range also 
shows the rapid variation of the individual parameters. Differential cross sections have 
been analyzed for the p-C 12 and p-O^ 6 system between 10 and 20 MeV (refs. 36 and 37). 
Cross-section and polarization data have been used in the analysis of p-N^ scattering 
at 20. 7 MeV (ref. 38) and p-O 16 scattering at 16. 3, 18. 5, and 21. 0 MeV (refs. 31 and 38). 

The rapid variation of the parameters may be due to the existence of intermediate 
structures which can yield broad resonances (a few hundred keV wide) in the cross sec- 
tion. Lowe and Watson (ref. 39) have taken these resonances into account by adding a 

contribution due to the coupling of three states in to the optical model scattering 

12 

amplitudes. Between 20 and 30 MeV, p-C cross sections and polarizations have been 
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analyzed by Tamura and Teresawa (ref. 40). The nonmonotonic variations of the cross 
section with energy could be reasonably well explained by the optical model if Breit- 
Wigner resonance terms were added. However, neither the polarizations nor the reac- 
tion cross sections are well represented by the parameters from the analysis of Tamura 
and Teresawa. 

For any further analysis of proton scattering by light nuclei in this "resonance 
region, ” it might be helpful to know how the optical potential would look if there were 
no resonances. 

An optical model analysis of the type described previously (covering a large range in 
nuclear masses) is one way to proceed. The "average" parameters of the optical poten- 
tial can be extrapolated into the region of low nuclear mass. 

The other possibility is an optical model analysis for a typical light nucleus over a 
large energy range. If at high enough energies a region can be found where the optical 

model is valid, those "average" parameters might be extrapolated to lower energies. 

12 

With the exception of C (which is a poor representative for light nuclei because of the 
strong coupling of the first excited state to the ground state), such a systematic study 
has not been done previously. 

For this reason Boschitz et al. (ref. 41) have extended the p-0 A polarization meas- 
urements at NASA (at 16. 3, 17. 8, and 21. 0 MeV) to 52. 5 MeV, using the 88-inch cyclo- 
tron at Berkeley. The polarization data have been supplemented by differential cross- 
section measurements at 42. 0 and 52. 5 MeV. The polarization data are shown in fig- 
ure 9. The lines drawn through the data points serve only to guide the eye. A signifi- 
cant change in the polarization pattern is found between 21 and 25 MeV. Above this 
energy, the polarization distributions vary slowly with energy. The analysis of these 
data has been performed in the following way: 

First, only those data were considered for which there appears to be a smooth 
energy dependence. The search was limited to seven parameters, which meant setting 
the geometrical parameters of the spin-orbit term equal to the corresponding para- 
meters of the central potential. The radius of the real central potential was succes- 
sively fixed at 1. 20, 1. 25, and 1. 30 F; the search code varied V, W, and V go , while 
simultaneously gridding on a, a^, and rj. For the search at 32. 8 MeV, the 30. 3 -MeV 
cross-section data from Harwell were used; while at 36. 7 MeV, interpolated cross sec- 
tions were used in the search routine. The best fits were selected by visual inspection 
under the constraint that the parameters were reasonable and compromising over the 
energy range. The best fits are seen in figure 10, and the corresponding parameters 
are tabulated in table IV. Clearly the agreement with the experimental data is not satis- 
factory, in particular, because all seven parameters are variables with energy. 

For this reason, radius r„ and diffuseness a o ^ of the spin-orbit potential and 
volume absorption W were added to the variables to be searched. The greatest im- 
provement in the fits to the 42. 0- and 52. 5 -MeV data was achieved by decreasing the 
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Polarization, 



Center of mass angle, 0 CM , deg 


Figure 9. - Energy dependence of polarization p-O 16 scattering. (Curves are faired through data points. Experimen- 
tal data taken from refs. 19, 20, and 41). 
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TABLE IV. - BEST CONSISTENT PARAMETERS WHEN A SEVEN- 
PARAMETER POTENTIAL IS USED 


Energy, 

Strength of 

Strength of 

Strength of 

Diffuseness 

Radius of 

Width of 

Radius of 

E, 

real poten- 

surface 

spin -orbit 

of real 

real 

surface 

surface 

MeV 

tial, 

part of 

potential. 

potential, 

potential, 

part of 

part of 


V 

potential, 

W D 

V so 

a = a 

so 

r = r 
o so 

potential, 

a i 

potential, 

r i 

52.5 

31.4 

3.6 • 

5.7 

0. 70 

1.20 

0.44 

1.25 

42.0 

36.0 

3.5 

8.1 

.75 

1.20 

.76 

1.30 

36.7 

41.3 

6.3 

8.1 

.67 

1.25 

. 50 

1.30 

32.8 

41.0 

5.4 

7. 5 

.65 

1.25 

. 50 

1.35 




Center of mass angle, 0 CM , deg 

Figure 11. - Best fits obtained to polarizations and differential cross sections for 0 16 with 10-parameter optical potential. 
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spin-orbit radius. Further improvements were obtained by decreasing a and adding 
volume absorption. The best fits obtained are shown in figure 11. By allowing the X 
values to exceed those of the best fits by a factor of about 2, quite systematic sets of pa- 
rameters were found (table V). 

These paremeters were extrapolated to the lower energies (25. 1, 23. 4, and 21. 1 
MeV), but there was no success in fitting the data. An extensive search showed that the 
parameters which do yield good fits become progressively inconsistent with decreasing 
energy. The best overall fits to polarizations and cross sections are shown in figure 12, 
and the corresponding parameters are given in table V. This table shows that in the 
energy range 30 to 50 MeV five geometrical parameters can be kept constant while the 
remaining five parameters vary smoothly with energy. The values of V and V__ de- 
crease almost linearly with energy. The energy dependence is approximately dV/dE ~ 
-1. 4 and dV g0 /dE » -0. 25. 

The total absorption (volume plus surface) stays nearly constant, but the volume ab- 
sorption increases with proton energy by the same amount as the surface absorption de- 
creases. A puzzling feature is the energy dependence of r . Many searches show 

SO 


TABLE V. - BEST CONSISTENT PARAMETERS WHEN A 10-PARAMETER POTENTIAL IS USED 


Proton 

energy, 

E, 

MeV 

Strength of 
real 

potential, 

V 

Strength of 
surface 
part of 
imaginary 
potential, 

W D 

Strength of 
volume 
part of 
imaginary 
potential, 
W 

Strength of 
spin- 
orbit 
potential, 

V so 

Diffuseness 
of real 
potential, 
a 

Radius of 
real 

potential, 

r o 

Width of 
surface 
part of 
potential, 

a i 

Radius of 
surface 
part of 
potential, 

r i 

Width of 
spin- 
orbit 
potential, 

a so 



Radius of 
spin- 
orbit 
potential, 

r so 

52. 5 

29.8 

3.9 

3. 0 

3.7 

0 . 

65 

1 . 

25 

0 . 

47 

1. 

30 

0.47 

1.05 

42. 0 

36. 1 

6.4 

1.0 

4.4 











1. 13 

36. 7 

41. 1 

7.2 

— 

7. 5 











1. 26 

32.8 

43.0 

7.8 

— 

8.0 











1. 22 

25. 1 

39.7 

4.7 

— 

4.3 










38 

1. 10 

23.4 

35.6 

4. 1 

— 

2.9 









. 33 

1 . 10 

21.0 

45.0 

10.7 

— 

2.7 






30 



.03 

1 . 08 


that it is impossible to obtain satisfactory fits over this energy range with a constant 
value of r go . At the lower energies the parameters are fairly inconclusive, especially 
in the case of a , which becomes unrealistically small at 21. 0 MeV. 

bO 

At the time when this analysis was made, no total reaction cross-section data were 

available in this energy range. Therefore, only a qualitative comparison could be made 

between theory and experiment by using total reaction cross-section data (refs. 42, 15, 
12 27 

and 43 to 46) for C and A1 (fig. 13). Recently the total reaction cross section has 
been measured for p-O 10 between 27 and 45 MeV (ref. 47). The predictions from this 
analysis are quite consistent w ith the data. 
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0 10 20 30 40 50 60 70 80 90 100 


Proton energy, Ep, MeV 

Figure 13. - Total reaction cross sections calculated with parameters in table V 
shown in comparison with experimental data for C 12 and A! 27 . 

The proton scattering from O ad is presently being studied between 20 and 30 MeV by 
Appel et al. (ref. 48). When the differential cross sections at large angles were meas- 
ured as a function of energy, resonances in the elastic scattering channel were observed 

at incident proton energies of 21. 4, 22. 3, 23. 2, and 27. 0 MeV. These resonances may 

17 

be related to compound states in F . 

COMPARISON WITH OPTICAL MODEL ANALYSES AT HIGHER ENERGIES 

For comparison, optical model analyses at higher energies will be briefly described. 
The discussion is restricted to analyses of polarizations and cross sections for several 
nuclei, so that average trends in the parameters can be well established. At 30 MeV, 
differential cross sections, polarizations, and total reaction cross sections have been 
measured (refs. 7 to 10) for a number of nuclei. Three analyses of these data have 
been reported. 

An optical model study using the differential cross-section data only was made by 
Barrett et al. (ref. 49) using the geometrical parameters previously determined by 
Perey in the analysis of 10- to 22 -MeV proton scattering data (i. e. , r Q = rj = r gQ = 

1. 25 F, a = a g0 = 0. 65 F, and a j = 0. 47 F). Barrett, et al. obtained satisfactory fits 
to the differential cross sections, except in the backward hemisphere. The polarizations 
predicted by these parameters, with the use of accepted values for the strength of the 
spin-orbit potential and the same geometrical parameters as for the real central poten- 
tial were in qualitative agreement with the experimental data. 

A more extensive analysis by Greenlees and Pyle (ref. 50) included polarizations as 
well as differential cross sections in the multiparameter search. A set of geometrical 
parameters (r Q = 1. 20 F, ^ = 1. 25 F, r gQ = 1. 10 F, and a = aj = a gQ = 0. 7 F) was found 
which gave good fits to both the differential cross section and the polarization data and 
simultaneously predicted the total reaction cross section correctly. The potential 
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strengths obtained were consistent with values of other energies. The larger value of the 
imaginary diffuseness (0. 7 F) in comparison with the values from the present analysis 
(0. 65 F) and Perey's work (0. 47 F) shows an increase of aj with energy. This corre- 
sponds to the increasing total reaction cross sections. Although Greenlees and Pyle 
found that good fits were obtained with a predominantly surface form for the imaginary 
central well, a slight improvement could be achieved by the addition of a small volume 
part. 

The most detailed analysis at this energy was done by Satchler (ref. 34). Satchler 
observed a smaller radius (r Q = 1. 12 F) and a larger diffuseness (a = 0. 75 F) for the 
real central potential than had previously been inferred. The symmetry term (N - Z)/A 
could be included either in the depth or in the radius of the real well. Satchler also sug- 
gested an (N - Z)/A term in the surface absorption potential. A volume absorption of 
about 3 MeV was required. Fits to the polarizations showed evidence for a spin-orbit 
coupling with a smaller radius and/or a smaller diffuseness than is used lor the real 
central potential. Satchler proposed two average potentials which gave an equally good 
overall account of the data. They are given in the following table: 


Potential 1 

V = 47. 5 + 0. 4 -5— + 30 --1^ 
AV3 A 

W n = 4. 5 + 16 ^-5 
U A 

r Q = 1. 12 F 


Potential 2 


V= 51 + 0.4-5— 
A 1 / 3 

W n = 4.25+ 16 
u A 

r n = 1. 09 + 0. 25 --^5 

U A 


W = 3 MeV; V SQ = 6. 1 MeV; i-j = 1. 33 F; r g0 = 1. 0 F; a = a g0 = 0. 75 F; 

fo. 58 F for Fe, Ni, Co, and Cu 
a^ = < 0. 65 F for Sn 
ip. 75 F for Pb 


Recently, at 40 MeV polarizations and cross sections have been measured (ref. 12) 
12 208 ' 
for 11 nuclei between C and Pb . The data were analyzed (refs. 12 and 35) in terms 

of the optical model, and the following " average parameters" were obtained: r = 

1. 16 F, rj = 1. 37 F, r go = 1. 06 F, a = 0. 75 F, ^ = 0. 63 F, a g0 = 0. 74 F, and °V go = 

6. 0 MeV. If the Coulomb correction term is taken to be the same as in Perey’s analy- 
sis, equal to 0.4 Z/A* 7 ^, the real central potential could be described by 
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N - Z 


V= 41. 1 + 0.4 -Z— + 26.4 

A 1 / 3 A 

The depth of the symmetry term in the potential, 26. 4 MeV, is in good agreement with 
Perey’s results (ref. 13). Volume, as well as surface absorption, is needed for good 
fits for nuclei heavier than calcium. The volume absorption varied between ~5 and 

7 MeV; the surface absorption increased with nuclear mass from W_, « 1. 0 MeV for 

54 208 u 

Fe to Wq « 4. 5 MeV for Pb . The average parameters inferred from this analy- 
sis also provide a fairly good fit to the 30-MeV data. By comparing the real well depth 
at both energies, an energy dependence for the real central well can be obtained; it is 
dV/dE = -0. 22±0. 03. This value is less than half of what was found between 10 and 
20 MeV. 


CONCLUSIONS 

An average potential was obtained which is able to predict quite well polarizations, 

differential cross sections, and total reaction cross sections for a large number of 

nuclei in the energy range from 16 to 22 MeV. The parameters of this potential are 

much closer to Perey's parameters than to the ones proposed by Rosen and coworkers. 

In particular, a dependence on the target mass is clearly established for the real central 

potential as well as for the depths and width of the surface peaked imaginary potential. 

For depth, diffuseness, and radius of the spin -orbit potential, the following average 

parameters were found: V ar , = 5. 5 MeV, a ari = 0. 55 F, and r ori = 1. 12 F. The smaller 

SO so so 

value of r in comparison with the radius of the real central well r = 1. 25 F is in 
so o 

good agreement with recent results for 30 and 40 MeV. There is no need for an imagi- 
nary part or a different form in the spin-orbit potential. 

Anomalously low values for the parameters of the spin-orbit potential were encoun- 
tered in the mass regions of doubly closed shells. No effect of this kind could be clearly 
established for the closure of only one shell. In the A = 40 region the anomaly is energy 

dependent and is most pronounced at energies where a new partial wave is entering. The 

40 40 

large difference in the polarization distribution for Ar and Ca near 17 MeV can be 

ascribed to a difference in the "effective strength, " the product of depth and diffuseness 

of the spin-orbit potential. A smooth energy dependence in this "effective strength" 

with some of the remaining parameters fixed (a = aj = 0. 65 F, r Q = r j = 1. 25 F, and 

r g0 = 1. 10 F) and some fluctuating slightly around average values (V^, V Ca « 47. 5 

MeV, W . ~ 7. 5 MeV, W„ « 5. 3 MeV) can predict the changing polarization and dif- 

Ar '-' a 40 39 40 

ferential cross section for Ar , K , and Ca between 14. 5 and 21. 0 MeV remarkably 
well. 
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The analysis of the energy dependence in p-O scattering shows that the elastic 

scattering of protons between 30 and 50 MeV is well described by the optical model. 

12 

This is in striking contrast to C where polarization and differential cross-section data 
have been analyzed in the comparable energy range by Craig and coworkers (29 and 
49 MeV) and by Fricke and coworkers (40 MeV). In none of these analyses was it pos- 
sible to obtain simultaneous fits to polarization and cross-section data. If one tries to 
extrapolate the optimum parameters for O between 30 and 50 MeV to lower energies, 
one meets with increasing difficulties in fitting the data. Fairly good fits can still be 
obtained at 25. 1 and 23. 4 MeV with a number of extrapolated parameters, but at 
21.0 MeV satisfactory fits can be obtained only with parameters which greatly differ 
from the extrapolated ones. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 8, 1968, 

129-02-04-06-22. 
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